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1.  INTRODUCTION 


In  weapon  testing  operations,  the  coordinates  of  the  burst  point  of  an  artillery  shell  are 
obtained  by  visual  observations  using  special  theodolites  that  are  placed  on  observation 
towers  (see  Roberts  1990  and  1991).  The  data  consist  of  observed  azimuths  and  elevations  of 
the  burst  point  and  the  ensuing  regression  problem  is  a  classic  problem  of  geodesy.  Roberts 
(1991)  describes  a  special  treatment  of  the  problem,  whereby  only  the  azimuth  observations 
are  used  to  compute  the  ground  coordinates  of  the  burst  point.  However,  the  elevation 
measurements  do  contribute  to  the  determination  of  the  ground  coordinates  and  should  not 
be  ignored.  This  report  describes  the  treatment  of  the  complete  data  set  in  accordance  with 
the  least-squares  principle  whereby  the  sum  of  the  squared  corrections  of  all  angle 
observations  is  minimized. 

2.  BURST-POINT  MEASUREMENT  PROCESS 

To  obtain  the  coordinates  of  the  burst  point  of  an  artillery  shell,  the  direction  to  the 
burst  is  measured  from  a  number  of  observation  towers  using  special  theodolites.  The 
observations  provide  the  azimuth  and  elevation  of  the  burst  point.  The  coordinates  of  the 
burst  point  can  be  computed  from  these  data  if  observations  from  at  least  two  towers  are 
available.  In  Section  3  we  present  a  least-squares  method  for  this  task.  The  result  of  the 
calculations  includes  estimates  of  the  burst-point  coordinates,  standard  deviations  of  the 
coordinates,  and  correlation  coefficients  between  the  coordinates.  A  typical  geometrical 
arrangement  of  observation  towers  and  weapon  is  such  that  sizable  correlations  between  the 
coordinate  estimates  can  be  expected.  Therefore,  it  is  important  to  have  estimates  of  the 
correlation  coefficients,  in  particular  if  an  average  burst  point  from  various  shots  is  to  be 
calculated,  because  the  variances  and  covariances  of  single-shot  coordinates  enter  into  the 
calculation  of  the  average  point. 

From  the  coordinates  of  the  burst  point,  coordinates  of  the  cannon,  and  the  firing 
direction  (defined  by  an  azimuth  angle  4>w),  we  compute  the  range,  deflection,  and  firing 
range  of  the  cannon.  These  quantities  are  defined  as  follows  (see  Figure  1): 

Range:  The  range  r  is  the  distance  between  the  cannon  and  the  projection  of  the 

burst  point  onto  the  level  plane.  The  projection  is  called  the  impact  point. 

Deflection:  The  angular  deflection  6  is  the  angle  between  the  firing  direction  and  the 
direction  from  the  cannon  to  the  impact  point.  The  metric  deflection  d  is 
defined  by 
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Impact  point 


Figure  1.  Definition  of  range  and  deflection 

d  =  r  sin£  .  (1) 

Firing  range:  The  firing  range  f  is  the  component  of  the  distance  (range)  to  the 
impact  point  in  the  firing  direction,  or 

/  =  r  cos  6  .  (2) 

3.  ESTIMATION  OF  BURST-POINT  COORDINATES 

Let  the  unknown  coordinates  of  the  burst  point  be  ( x,y,z )  and  the  coordinates  of  the 
<-th  observation  tower  be  (xt,yt,zt).  Let  the  observed  azimuth  and  elevation  be  <j>t  and  $t, 
respectively,  and  let  T  be  the  number  of  observation  towers.  Then  the  following  relations 
can  be  read  from  Figure  2: 


Figure  2.  Tower  and  burst-point  coordinates. 

9t(^t>x>v)  —  (*  “  xt)  tan^t  —  (y  —  y<) =  0  ,  t  =  1,2,..,  T  ,  (3) 

and 

A|(^;*>y^)as(^-^)2-[(*-*»)2  +  (y-y<)2]tan2fl<=o  ,  *  =  i,2,..,r.  (4) 

Equations  (3)  and  (4)  are  the  model  equations  of  the  regression  problem.  The  regressands  are 
the  angles  <f>t  and  0t,  t  =  1  ,...,T,  the  regressor  variables  are  the  tower  coordinates  (xt,yt,zt), 
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and  the  free  model  parameters  are  the  burst-point  coordinates  ( x,y,z ).  Let  P,  be  estimated 
variance-covariance  matrices  of  the  observations,  and  let  c+t  and  ctt  be  the  corrections 
(residuals)  of  the  observations  j>t  and  6t,  respectively.  Then  the  least-squares  regression 
problem  can  be  formulated  as  follows: 

Minimize  W  =  E  (c^t,ctt)  Pj-1  °*1  ,  (5) 

subject  to  gt(4t  +  c^;  x,y)  ~  0  , 

and  ht(0t  +  c,t\  x,y,z)=0  , 


This  is  a  least-squares  problem  with  constraints  in  the  form  of  simultaneous  equations.  It 
can  be  solved,  for  instance,  with  the  help  of  the  utility  routine  COLSMU  (Celmi^S  1979). 
However,  in  the  present  case,  some  simplifications  apply  that  allow  one  to  use  the  simpler 
routine  COLSAC  (Celmipg,  l.c.)  for  least-squares  problems  with  scalar  constraints.  First,  we 
can  assume  that  all  angle  observations  are  independent  and  have  equal  accuracies. 
Therefore,  the  variance-covariance  matrices  Pt  in  Eq.  (5)  can  be  set  equal  to  unit  matrices 
and  the  objective  function  defined  by 

W-W,+  Wt-  £  cl  +  £  cl  .  (7) 

Second,  because  every  equation  in  (6)  contains  only  one  scalar  observable,  the  constraints  can 
be  defined  as  the  following  2  T  scalar  equations: 

M  +  c*i!  *,y)=0  ,  <=1,2,..., T  ,  (8a) 

At_T(0|_r  +  c,|_T;  r)=0,  t  =  T+l,T+2 . 2 T.  (8b) 

Eqs.  (7),  (8a),  and  (8b)  define  a  least-squares  problem  with  2  T  observations,  2  T  scalar 
constraints,  and  three  parameters. 

Roberts  (1991)  suggests  for  the  solution  of  the  regression  problem  a  method  that  is 
different  from  the  outlined  approach  and  does  not  provide  the  least-squares  solution.  To 
show  the  difference  between  a  least-squares  solution  and  Roberts’  method,  we  first  derive  a 
set  of  normal  equations  for  the  least-squares  problem.  We  start  by  solving  the  constraint 
equations  (6)  for  the  unknown  residuals.  The  results  are 


c  arctan 


y-Vt 

X  —X( 


cet  —  ~  a-rctan 


_ *-*t _ 

[  (x  —  xt)2  +  (y  —  y<)2  J1/2 


<=1,2,...,  T  , 

<=1,2,..., T  . 


(9) 
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(10) 

(11) 


Now  the  only  unknowns  in  the  objective  function  W  =  W1  +  W2  are  the  free  model 
parameters  x,  y,  and  z.  We  obtain  equations  for  these  parameters  by  setting  equal  to  zero 
the  partial  derivatives  of  W  with  respect  to  the  unknowns.  The  ensuing  system  of  normal 
equations  is 

dW  _  dWx{x,y)  dW^x,y,z) 

dx  dx  dx 

dW  _  dWfay)  dW£x,y,z) 

dy  dy  dy 

dW  _  aiy2(z,y,z) 

dz  dz 

The  least-squares  values  of  x,  y,  and  z  are  solutions  of  this  system  of  equations.  Numerically 
solving  Eq.  (12)  is  equivalent  to  solving  the  least-squares  problem  (7)  and  (8)  with  the  utility 
program  COLSAC. 

Roberts  (1991)  solves,  instead  of  the  normal  equation  system  (12),  the  following  simpler 
system: 

dWx(x,y) 
dx 

dW\{x,y) 

dy 

dW2(x,y,z) 
dz 

This  modification  of  the  normal  equation  system  is  justified  if  one  assumes  that  the 
dependence  of  the  azimuth  angle  0  on  the  altitude  z  is  negligible.  The  assumption  might  be 
true  for  typical  firing  tests,  but  sufficient  conditions  for  the  assumption  to  hold  have  not  been 
elaborated.  If  the  assumption  is  not  true,  then  the  solution  of  the  modified  set,  Eq.  (13),  is 
different  from  the  least-squares  solution.  The  magnitude  of  the  difference  is  not  known 
Hence,  because  the  complete  set  of  normal  equations,  Eq.  (12),  provides  the  least-squares 
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solution  under  all  conditions,  we  see  no  reason  to  use  the  modified  set.  Roberts  does  not  give 
a  reason  for  the  neglect  of  the  terms  dW2/dx  and  dW2/dy.  Presumably,  the  terms  were 
removed  to  facilitate  the  numerical  solution.  However,  the  first  two  equations  of  the 
system  (13)  constitute  a  system  of  simultaneous  non-linear  equations  for  x  and  y  that  cannot 
be  formally  solved.  Roberts  obtained  a  numerical  solution  of  that  system  by  linearization 
and  iteration.  He  could  have  used  the  same  numerical  technique  on  the  complete  equation 
system  (12)  and  obtained  the  correct  least-squares  solution. 


In  the  example  presented  in  Section  5,  the  numerical  solution  of  the  modified  Eq.  (13)  is 
found  to  be  different  from  the  numerical  solution  of  the  complete  Eq.  (12). 


A  numerical  solution  of  the  general  regression  problem  can  be  obtained  also  with  the 
help  of  commercial  software  for  non-linear  data  fitting.  Such  a  software  typically  requires 
that  each  constraint  equation  contains  exactly  one  regressand  and  that  the  constraints  are 
explicitly  solved  for  the  observations;  that  is,  they  must  be  of  the  type 

*<-/<«.*)  +  <*.  (M) 

where  are  the  observations  (the  "regressand  variables"),  X{  are  fixed  parameter  vectors 
(the  "regressor  variables"),  f)  is  a  free  parameter  vector  of  the  regression  functions  Ji ,  and 
c^i  are  the  corrections  (residuals)  of  the  observations.  In  the  present  case,  this  form  of 
constraints  can  be  obtained  by  solving  Eqs.  (9)  for  the  observations.  (The  utility  routine 
COLSAC  does  not  mandate  any  particular  constraint  fc  mutation  and  accepts  constraints  in 
the  general  form  /,•  (^,+e^, ,  X{ ,  0)  —  0.  Therefore,  we  used  the  simpler  implicit 
formulations  (3)  and  (4)  instead  of  Eqs.  (9)  in  our  numerical  calculations  with  COLSAC.  )  As 
an  example  of  a  commercial  software,  we  chose  the  program  3R,  Release  1990,  of  BMDP 
Statistical  Software,  Inc.  To  comply  with  3R,  we  expressed  the  constraints  in  the  form  of  a 
regression  model  as  follows: 


4>t  »|  +  6t  i2  — 

[y  —  yt  -I  r 
arctan - J  *i  +  L  Cretan 


z  —  z. 


X  —  X, 

+  Cft  t'j  +  ctt  »2 


[  (*  -  *<)2  +(y  -  vt?  ]1/2 
t  -  1,  2,  .  .  .  ,  T. 


]  *2  + 


(15) 


The  quantities  and  i2  are  called  indicator  variables.  Their  values  are  *j  =  1  and  i2  =  0  for 
azimuth  measurements  and  *j  =  0  and  i2  =  1  for  elevation  measurements.  The  regressor 
variables  are  the  tower  coordinates  yt,  and  zt  ,  and  the  components  of  the  free  model 
parameter  vector  are  the  burst-point  coordinates  x,  y,  and  z. 


-5  - 


4.  COMPUTATION  OF  RANGE  AND  DEFLECTION 


Let  ( x,y,z )  be  the  coordinates  of  the  burst  point,  {xw,yw,zw)  be  the  coordinates  of  the 
gun,  <j>w  be  the  azimuth  angle  of  the  firing  direction,  and  7?  be  a  unit  vector  in  the  direction 
of  firing  defined  by 


’  , 

COB<f>w 

sin^, 


(16) 


We  define  a  range  vector  by 


r  = 


x-x. 


y-yy> 

The  range  is  the  length  of  the  range  vector,  that  is, 


The  metric  deflection  is 

d  =  [  nX~r  ],  =  -(*-  xm)  sin^  +  [y  -  yw)  cos <f>w  . 
The  angular  deflection  is 


6  =  arcsin(d/r) 


and  the  firing  range  is 

/  =  ft  ■  7  —  r  cos  (5  =  (x  -  xw)  cos<j>K  +  (y  —  yw)  sin  <f>v  . 


(17) 


(18) 

(19) 

(20) 
(21) 


To  obtain  accuracy  estimates  of  the  ranges  and  deflections,  we  use  the  variance 
propagation  formula.  Let  PB  be  the  variance-covariance  matrix  of  the  burst-point 
coordinates.  It  is  a  3X3  matrix  and  one  of  the  outputs  of  the  COLSAC  utility  routine.  Then 
the  variance-covariance  matrix  of  the  range,  metric  deflection,  and  height  of  burst  is 


Pm  = 


d(r,d,z) 

d{x,y,z) 


Pb 


[  9[r,d,z) 
[  d{x,y,z) 


T 


(22) 


Computing  the  derivatives  one  obtains  the  explicit  formula 


{x  -  xj/r  (y  -  yw)/r  0 

(x  -  xj/r 

—  sin  <f>w 

0 

p*- 

-sin^  cos^w  0 

Pb 

(y  -  yJA 

cos  <t>v 

0 

0  0  1 

0 

0 

r2 

(23) 


The  variance-covariance  matrix  of  the  range,  angular  deflection,  and  height  of  burst  can  be 
obtained  in  the  same  manner  and  is 
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(24) 


(x  -  x  J/r 

■S' 

1 

o 

(x  -  xj/r  -  (y  -  yj/r2  0 

-(y-y„)A2 

O 

N 

W 

1 

H 

Pb 

(y  -  yJA  (*  -  *JA2  o 

0 

0  1 J 

o 

o 

The  formula  for  the  variance-covariance  matrix  of  the  firing  range,  metric  deflection,  and 
height  of  burst  is 


[  C06^_ 

sin 

0 

COS 

—  sMw 

0 

» 

.a 

to 

1 

cos  <f>w 

0 

Pb 

sin  4>w 

cos^„ 

0 

0 

0 

1 

0 

0 

1. 

(25) 


The  formula  for  the  variance-covariance  matrix  of  the  firing  range,  angular  deflection,  and 
height  of  burst  is 


cos  <f>w 

sin^w  0 

C<xs4>u 

-  (y  -  yJA2  o 

(y  yJA2 

(x  -  xjA2  0 

Pb 

sin  4>m 

(x-xj/r2  0 

0 

0  1  , 

0 

0  1 

(26) 


5.  EXAMPLES 

We  use  the  examples  from  Roberts  (1991).  The  coordinates  of  four  observation  towers 
and  a  gun  are  given  in  Table  1.  (We  have  subtracted  2000  m  from  the  x-coordinates  and 
10000  m  from  the  y-coordinates  for  simplicity.  Also,  the  azimuth  of  the  firing  direction  is 
modified  by  subtracting  35°  because,  according  to  Roberts  (1991),  a  special  coordinate  system 
was  used  for  the  measurements.)  Table  2  lists  observations  of  azimuths  and  elevations  for 
eight  rounds  labeled  11  through  18.  The  azimuths  are  again  reduced  by  35°.  Results  of 
regression  by  the  least-squares  method  (Eqs.  (7)  and  (8),  or  the  equivalent  Eqs.  (10),  (11),  and 
(12) )  that  were  obtained  with  the  utility  program  COLSAC  are  shown  in  Tables  3  and  4. 
Table  3  contains  the  estimated  burst-point  coordinates,  translated  to  a  coordinate  system 
with  the  cannon  location  (x„,yw,z„)  at  the  origin,  and  estimates  cv  ef,  and  e,  of  the 
standard  deviations  of  the  burst-point  coordinates.  All  entries  in  Table  3  are  in  metres. 
Table  4  contains  estimated  correlation  coefficients  between  the  burst-point  coordinates. 

The  variance-covariance  matrix  Pb  of  the  coordinates  is  computed  from  the  standard 
deviations  and  the  correlation  coefficients  Cg  by  the  formula 

Pik  =  ei  ek  c*  ,  (27) 

where  p #  is  an  element  of  the  matrix  PB  and  the  indices  t  and  k  may  take  the  values  x,  y,  or 
2. 


Identical  numerical  results  were  obtained  with  the  utility  program  3R,  Release  1990,  of 
BMDP  Statistical  Software,  Inc.,  that  uses  the  constraints  (15).  This  software  also  provides 
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tests  for  the  dependence  of  the  parameters  x,  y,  and  2  on  data  deviations.  The  tests 
indicated  that  the  dependence  is  essentially  linear  for  the  ranges  of  interest.  This  finding 
justifies  the  use  of  the  linearized  law  of  variance  propagation  for  the  parameter-variance 
estimates  in  the  utility  routine  COLSAC  and  for  the  computation  of  deviations  of  range  and 
deflection  by  the  formulas  of  Section  4. 

For  comparison,  Tables  5  and  6  display  results  obtained  by  solving  the  equation 
system  (10),  (11),  and  (13)  that  correspond  to  Roberts’  solution  and  agree  with  the  values 
reported  by  Roberts  (1991).  We  have  supplemented  his  solution  with  estimates  of  coordinate 
standard  deviations  and  correlation  coefficients. 

The  estimates  of  coordinates  by  the  two  methods  differ  only  by  less  than  3  m, 
indicating  good  quality  and  consistency  of  the  data.  (If  the  data  contain  large  measurement 
inaccuracies,  then  the  differences  between  results  from  different  estimation  methods  can  be 
significant.)  However,  the  estimated  standard  deviations  vary  significantly.  For  instance, 
consider  the  estimates  of  the  standard  deviations  of  angle  observations  that  are  listed  in 
Table  7.  From  the  measurement  technique,  we  expect  equal  uncertainties  of  azimuth  and 
elevation  angle  measurements,  respectively.  Consequently,  in  our  analysis  we  assumed  equal 
accuracies  and  the  analysis  produced  for  each  round  one  estimate  of  the  standard  deviation 
of  all  angle  measurements.  In  contrast,  Roberts  obtained  different  error  estimates  for 
azimuth  and  elevation  readings,  respectively,  because  he  used  Eq.  (13)  instead  of  the  least- 
squares  normal  equations  (12).  The  consequences  are  most  apparent  for  Rounds  11,  15,  and 
16  with  larger  data  scatter.  It  seems  that  Roberts’  method  arbitrarily  assigns  large  errors 
either  to  azimuth  or  to  elevation  observations  making  the  other  observations  appear 
extremely  accurate.  We  fail  to  see  any  justification  for  such  an  assignment  of  residuals. 

Tables  8  through  11  list  the  results  in  terms  of  firing  range,  metric  deflection,  and 
height  of  burst.  The  quantities  in  the  table  are  given  in  metres  and  they  were  computed 
from  the  coordinates  of  the  burst  points  as  described  in  Section  4.  The  results  are  also 
shown  in  Figures  3  through  8. 

Figures  3,  4,  and  5  show  the  locations  of  all  burst  points  calculated  by  the  least-squares 
method.  Figures  6,  7  and  8  show  the  corresponding  results  obtained  by  Roberts’  method. 
The  accuracies  of  the  burst-point  coordinates  are  indicated  by  corresponding  one-standard- 
error  ellipsoids.  We  observe  in  Figure  3  that  Rounds  11  and  16  (denoted  by  "1"  and  "6", 
respectively)  stand  out  as  less  accurate.  In  contrast,  Figure  6  shows  that  Roberts’  method 
assigns  only  to  Round  11  large  errors  of  ground  coordinates,  whereas  the  location  error 
estimates  of  all  other  burst  points  are  overly  optimistic. 
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Figure  3.  Firing  range  and  metric  deflection 
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Figure  4.  Firing  range  and  height 

6.  TESTS  FOR  OUTLIERS 

Ad  important  part  of  the  treatment  of  burst-point  observations  is  the  detection  of 
gross  observational  errors.  A  theoretical  basis  for  the  detection  with  a  prescribed  confidence 
level  is  not  available  because  the  problem  is  not  linear  and  outlier  detection  theories  have 
been  elaborated  only  for  linear  regression.  One  may  therefore  use  statistical  tests  based  on 
linear  regression  as  approximations  or  establish  ad  hoc  tests  specifically  for  the  problem  at 
hand.  We  shall  discuss  both  approaches. 

Table  12  lists  the  standard  angular  deviations  e^  and  the  largest  residuals  of  the  angie 
observations  for  each  round  (e^  denotes  the  estimate  of  the  standard  angular  deviation  for  <f> 
and  0,  and  denotes  the  residual  of  an  angle  observation).  It  is  obvious  from  this  list  that 
Rounds  11  and  16  stand  out  as  candidates  with  bad  observations  (see  also  Figure  3). 
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Figure  5.  Metric  deflection  end  height 


Figure  0.  Firing  range  and  metric  deflection  by  Roberta 

Another  leas  pronounced  outlier  might  be  among  the  observations  of  Round  15.  A  perusal  of 
Table  12  shows  that  good  indicators  for  rounds  with  excessive  residuals  are  the  magnitudes 
of  and  |c^|max.  One  can,  for  instance,  postulate  for  the  present  arrangement  of  towers 
and  cannon  the  following  test  for  a  set  with  outliers: 

e*  >  0.20°  or  |c*lm„  >  0-25°  •  (28) 

For  an  outlier  itself,  one  might  use  the  same  criterion: 

|c*|>0.25°.  (29) 

To  establish  tests  of  this  type,  one  must  know  what  sizes  of  deviations  to  expect  in  a 
normal  operation,  that  is,  the  characteristics  of  the  distribution  of  the  observation  errors. 
Thus,  the  conditions  (28)  and  (29)  were  derived  by  comparing  results  from  the  eight  different 
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Figure  7.  Firing  range  and  height  by  Roberta 
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Figure  8.  Metric  deflection  and  height  by  Roberts 

rounds  in  Table  12.  Without  the  possibility  to  compare  between  rounds  (if  only  one  or  two 
rounds  have  been  measured  with  the  same  configuration  of  towers  and  cannon)  we  may  use 
the  third  column  of  Table  12  as  test  statistic.  It  is  obvious  that  this  statistic  requires  a 
careful  fine  tuning  of  the  outlier  test,  because  a  large  error  in  one  measurement  increases  the 
residual  as  well  as  e^,  and  makes  their  ratio  less  distinct  from  those  of  error-free 

measurements.  After  inspecting  the  data  in  the  third  column,  one  might  postulate  in  the 
present  case  the  following  condition  as  an  error  indicator: 

Mn»x/«*>  1-70  .  (30) 

Linear  regression  theory  derives  a  similar  test.  For  a  5%  confidence  level  with  eight 
observations  and  three  parameters,  an  outlier  is  indicated  if  (see  Barnett  and  Lewis  1978, 
p.  262  ff.) 
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(31) 


I  c+\  max  /  >  2.10  . 

Note  that  according  to  this  linear  theory,  all  observations  would  be  accepted. 

Outlier  detection  can  be  also  based  on  repeated  regression  with  one  observation  deleted 
in  turn  (see  Barnett  and  Lewis  1978,  p.  244  if.).  This  approach  is  not  feasible  if  the  number 
of  observations  is  large,  but  with  the  present  eight  observations  per  round,  only  a  moderate 
amount  of  additional  calculations  for  the  test  are  needed.  The  advantage  of  this  test  is  that 
it  makes  a  single  outlier  stand  out  very  pronounced.  For  illustrative  purposes,  we  first  show 
in  Table  13  all  relative  residuals  of  the  Rounds  11,  15,  and  16,  that  is,  of  those  rounds  that 
are  suspect  of  containing  outliers  according  to  the  test  (30).  The  observations  Round  11- 
Tower  8-Azimuth  and  Round  16- Tower  2-Elevation  clearly  stand  out  as  suspect.  The 
observation  Round  15-  Tower  8-Elevation  is  a  borderline  case  and  it  should  be  kept  in  the 
data  set  unless  there  are  other  indications  of  errors.  We  now  repeat  the  regressions,  deleting 
one  observation  in  turn.  The  idea  behind  this  approach  is  that  the  estimated  standard 
deviation  is  significantly  smaller  if  an  outlier  is  removed  from  the  set  but  changes  only  little 
if  a  good  observation  is  removed.  It  can  be  shown  that  in  linear  problems  this  approach  is 
equivalent  to  the  maximum  residual  statistics  test,  Eq.  (31),  (see  Barnett  and  Lewis  1978, 
p.  265),  but  it  accentuates  the  dichotomy  between  good  and  bad  observations.  Table  14  lists 
the  test  statistics  of  the  regressions  by  deleting  one  observation.  Each  entry  in  the  table  is 
the  ratio  of  the  standard  deviation  of  the  full  set  to  the  standard  deviation  e+  that  is 
obtained  if  the  corresponding  observation  is  removed  from  the  set.  Hence,  the  entries  are  the 
relative  magnitudes  of  the  standard  deviation  estimates  if  the  corresponding  observations  are 
left  in  the  data  set.  A  large  entry  indicates  a  possible  outlier.  The  suspect  observations  of 
Round  11  and  16  with  standard  deviation  ratios  of  5.3  and  12.1,  respectively,  stand  out  more 
prominently  in  Table  14  than  in  Table  13.  Outliers  might  be  identified  in  the  present 
regression  problem,  for  instance,  by  the  ad  hoc  test 

y*  /  <!*  >  3.  (32) 

For  Round  15,  there  is  no  prominent  deviation  from  unity  in  Table  14.  One  might  check  the 
observation  Round  15-Tower  8-Elevation  for  correctness,  but  one  would  not  discard  the 
measurement  based  only  on  the  Table  14.  We  have  no  such  hesitations  for  the  removal  of 
the  two  outliers  for  Rounds  11  and  16. 

After  removing  the  outliers  from  the  observations  of  Rounds  11  and  16,  we  obtain  new 
burst-point  coordinates  and  deviation  estimates.  The  new  values  are  listed  in  Tables  15  and 
16  and  shown  in  Figures  9,  10,  and  11.  The  largest  changes  are  in  the  deflection  of  Round  11 
that  has  changed  from  —1.42  m  to  -fll.39  m  and  in  the  height  of  Round  16  that  has 
decreased  by  14.17  m.  The  overall  appearance  of  the  burst-point  dispersion  has  improved  as 
indicated  by  a  comparison  of  Figure  3  with  Figure  9. 
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Figure  0.  New  firing  range  and  metric  deflection 


Firing  range,  m 

Figure  10.  New  firing  range  and  height 

7.  SUMMARY  AND  CONCLUSIONS 

Measurements  for  the  determination  of  burst-point  coordinates  consist  of  azimuth  and 
elevation  angles.  We  compute  the  coordinates  by  minimizing  the  sum  of  the  squares  of  the 
residuals  of  all  angle  observations.  The  result  is  a  set  of  burst-point  coordinates  that  agrees 
best  with  the  observations  in  a  least-squares  sense.  In  this  report,  we  describe  the  formalism 
for  the  calculation  of  the  coordinates  and  propose  two  utility  routines  for  the  numerical 
solution  of  the  regression  problem.  We  also  discuss  methods  for  the  detection  of  outliers  in 
the  data  sets. 

The  described  method  for  the  determination  of  burst-point  coordinates  differs  from 
that  of  Roberts  (1991).  His  method  does  not  minimize  the  sum  of  all  squared  residuals  and 
therefore  is  not  a  least-squares  algorithm  in  the  usual  sense.  We  show  and  discuss  in  an 
example  the  differing  results  that  are  obtained  by  the  two  algorithms.  We  conclude  that 
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Figure  11.  New  metric  deflection  end  height 

Roberts’  algorithm  should  not  be  used  except  when  accurate  burst-point  coordinates  are  not 
needed. 

In  discussing  tests  for  outliers  we  observe  that  teste  that  are  based  on  the  absolute  size 
of  residuals  require  the  least  additional  computing.  However,  they  have  the  drawback  that 
the  threshold  for  outlier  detection  depends  on  the  setup  of  the  experiments.  For  instance,  a 
radical  change  of  the  firing  range  or  of  the  arrangement  of  observation  towers  would  likely 
require  a  change  of  the  constants  in  the  tests  (28)  and  (30).  A  test  that  is  less  sensitive  to 
changes  in  the  experiment  is  based  on  repeated  regressions  by  deleting  one  observation  in 
turn.  It  requires  more  computation,  but  the  increase  of  computing  time  is  moderate  as  long 
as  the  number  of  observation  towers  is  not  excessive.  We  propose  to  develop  a  dedicated 
utility  program  that  includes  this  test  for  the  computation  of  burst-point  coordinates.  Such 
a  program  could  provide  coordinate  estimates  and  outlier  tests  likely  in  real  time. 
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Table  1.  Tower  and  cannon  coordinates. 


z,  m 

y,  m 

z,  m 

Tower  1 

1657.607 

2801.626 

ISBW| 

Tower  2 

1868.530 

3870.487 

Tower  3 

1754.335 

3094.383 

4.77 

Tower  4 

883.790 

1998.310 

12.13 

Gannon 

1812.273 

isa 

Firing  direction  azimuth  4>w  —  6° 

26' 
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Table  2.  Observed  azimuths  and  elevations. 
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Table  3.  Estimates  of  burst-point  coordinates. 


Round 

x-xw 

■S9 

mm 

z~zw 

e* 

m 

2530.75 

42.65 

283.92 

10.62 

952.65 

17.43 

— 

2408.26 

11.27 

293.79 

2.90 

852.55 

4.44 

R9. 

2547.62 

10.84 

296.12 

2.66 

952.83 

4.41 

14 

2388.24 

5.27 

292.72 

1.37 

839.60 

2.07 

15 

2529.40 

16.39 

286.87 

4.06 

303.23 

3.90 

16 

2570.77 

49.19 

292.70 

11.97 

293.64 

11.48 

17 

2501.56 

8.87 

284.22 

2.22 

297.19 

2.12 

18 

2533.77 

6.80 

292.17 

1.67 

274.36 

1.58 

Table  4.  Correlation  coefficients  to  Table  3. 


Round 

c*v 

c» 

c. . 

mm 

-0.45766 

0.80939 

-0.39097 

'MSm. 

-0.45865 

0.78380 

-0.38057 

m m 

-0.44965 

0.81068 

-0.38446 

14 

-0.45897 

0.77999 

-0  37980 

15 

-0.46651 

0.42803 

-0.21035 

16 

-0.46006 

0.41707 

-0.20001 

17 

-0.46872 

0.42178 

-0.20888 

18 

-0.46006 

0.39394 

-0.19135 

Table  5.  Roberts’  estimates  of  burst-point  coordinates. 


Round 

x-xw 

e. 

■5S 

er 

z—zw 

mm 

2532.65 

67.28 

283.58 

16.63 

953.28 

22.77 

■p* 

2410.52 

10.83 

293.26 

2.77 

853.26 

6.62 

2544.35 

7.86 

296.54 

1.92 

951.75 

6.51 

14 

2387.91 

3.97 

292.89 

1.03 

839.50 

3.18 

15 

2530.91 

2.70 

286.61 

0.67 

303.38 

9.00 

16 

2570.81 

2.59 

292.29 

0.63 

293.65 

26.75 

17 

2501.98 

3.93 

284.16 

0.98 

297.24 

4.75 

18 

2533.31 

7.90 

292.22 

1.94 

274.32 

2.64 
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Table  8.  Correlation  coefficients  to  Table  4. 


Round 

c*v 

c»* 

mm 

-0.46500 

0.97787 

-0.47923 

-0.46531 

0.50555 

-0.24874 

Hi 

-0.45680 

0.39813 

-0.19151 

14 

-0.46545 

0.38257 

-0.18859 

15 

-0.46741 

0.03052 

-0.01504 

16 

-0.46074 

0.00942 

-0.00455 

17 

-0.46958 

0.08326 

-0.04131 

18 

-0.46079 

0.27395 

-0.13326 

Table  7.  Standard  deviations  of  angle  observations. 


Least  Sq. 

Roberts 

Round 

and  Cg 

U 

22.34' 

34.76' 

5.18' 

6.37' 

6.04' 

6.63' 

wM 

5.61' 

4.02' 

6.54' 

14 

3.02' 

2.52' 

3.45' 

15 

8.48' 

1.39' 

10.90' 

16 

24.79' 

1.30' 

32.00' 

17 

4.68' 

2.07' 

5.80' 

18 

3.51' 

4.07' 

3.08' 

Table  8.  Estimates  of  firing  range,  metric  deflection,  and  height  of  burst 


Round 

Firing  range  / 

ef 

Deflection  d 

U 

Height  z 

ra 

2546.63 

41.85 

-1.42 

13.43 

952.65 

la 

2426.02 

11.05 

22.10 

3.64 

852.55 

2564.76 

10.64 

8.80 

3.37 

952.83 

14 

2406.01 

5.17 

23.28 

1.71 

839.60 

15 

2545.61 

16.08 

1.65 

5.15 

303.23 

3.90 

16 

2587.38 

48.28 

2.81 

15.24 

293.64 

11.48 

17 

2517.65 

8.70 

2.14 

2.81 

297.19 

2.12 

18 

2550.54 

6.67 

6.43 

2.12 

274.36 

1.58 
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Table  9.  Correlation  coefficients  to  Table  7 


Round 

cfd 

eU 

CiM 

mm 

-0.69764 

0.80855 

-0.59528 

im 

-0.69136 

0.78291 

-0.57323 

mm 

-0.69529 

0.80987 

-0.59373 

14 

-0.69018 

0.77907 

-0.57033 

15 

-0.70427 

0.42759 

-0.31729 

16 

-0.70357 

0.41672 

-0.30701 

17 

-0.70345 

0.42133 

-0.31204 

18 

-0.70133 

0.39353 

-0.29112 

Table  10.  Roberts'  estimates  of  firing  range,  metric  deflection,  and  height  of  burst. 


Round 

Firing  range  / 

ef 

Deflection  d 

e4 

Height  z 

mm 

2548.47 

66.01 

-1.98 

21.12 

953.28 

22.77 

n  -- 

2428.20 

10.62 

21.32 

3.49 

853.26 

6.42 

wm 

2561.55 

7.72 

9.59 

2.44 

951.75 

6.51 

14 

2405.69 

3.89 

23.48 

1.29 

839.50 

3.18 

15 

2547.08 

2.64 

1.23 

0.85 

303.38 

9.00 

16 

2587.37 

2.54 

2.40 

0.80 

293.65 

26.75 

17 

2518.06 

3.85 

2.03 

1.24 

297.24 

4.75 

18 

2550.10 

7.75 

6.53 

2.47 

274.32 

2.64 

Table  11.  Correlation  coefficients  to  Table  10. 


Round 

cji 

CU 

ci. 

n 

-0.70343 

0.97687 

-0.72421 

El 

-0.69677 

0.50499 

-0.37219 

HE 

-0.70079 

0.39774 

-0.29351 

14 

-0.69525 

0.38213 

-0.28140 

15 

-0.70503 

0.03049 

-0.02266 

16 

-0.70430 

0.00941 

-0.00696 

17 

-0.70413 

0.08317 

-0.06183 

18 

-0.70188 

0.27367 

-0.20259 
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Table  12.  Largest  residuals. 


Round 

e*.  degr. 

lc*L«.  de«r- 

1  c^l  m*x/et 

mm 

1 

0.684 

MEM 

If 

Bai'Sl 

0.154 

msm. 

0.093 

0.136 

1.459 

14 

0.050 

0.068 

1.351 

15 

0.141 

0.223 

1.577 

16 

0.413 

0.786 

1.903 

17 

0.078 

0.120 

1.540 

18 

0.058 

0.072 

1.226 

Table  IS.  Relative  residuals  |  |  /  . 


Round  11.  e^=  0.772° 

Tower  1 

Tower  2 

Tower  3 

Tower  4 

Azimuth 

Elevation 

-0.6996 

0.1848 

-0.8267 

-0.0391 

1.8365 

-0.2917 

-0.5456 

0.1897 

Round  15.  0.141° 

Tower  1 

Tower  2 

Tower  3 

Tower  4 

Azimuth 

Elevation 

-0.1041 

1.2984 

0.0282 

-0.3067 

-0.0402 

-1.5768 

■ 

Round  16.  =  0.413° 

Tower  1 

Tower  2 

Tower  3 

Tower  4 

Azimuth 

Elevation 

0.0496 

0.7996 

0.0436 

-1.9027 

-0.0406 

0.7433 

HI 
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Table  14.  Ratios  of  standard  deviations  c+  f  . 


Round  11. 

Tower  1 

Tower  2 

Tower  3 

Tower  4 

Azimuth 

Elevation 

0.971 

0.899 

1.477 

0.895 

5.323 

0.905 

■IfB 

givnrTTH 

Round  15. 

Tower  1 

Tower  2 

Tower  3 

Tower  4 

Azimuth 

Elevation 

0.896 

1.224 

0.895 

0.906 

0.895 

1.641 

0.904 

0.976 

Round  16. 

Tower  1 

Tower  2 

Tower  3 

Tower  4 

Azimuth 

Elevation 

0.895 

0.986 

0.895 

12.143 

1  '•  _ 

■ 

Table  16.  New  estimates  of  burst-point  locations. 


Round 

Firing  range  / 

«/ 

Deflection  d 

u 

Height  z 

11  old 

2546.63 

41.85 

mgm 

13.43 

952.65 

17.43 

22.34' 

11  new 

2552.04 

7.92 

— 

2.74 

953.38 

3.29 

4.20' 

16  old 

2587.38 

48.28 

1!M 

15.24 

293.64 

11.48 

24.79' 

16  new 

2587.60 

3.98 

1.26 

279.47 

1.07 

2.04' 

Teble  16.  New  correlation  coefficients  to  Table  16. 


j  Re  and 

cn 

CJ * 

ciz 

SB 

-0.61295 

0.80925 

-0.53645 

wM 

-0.70370 

0.34778 

-0.25203 
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